on the theory of supersolidity 
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We present a microscopic, many-body argument for supersolidity. We also illustrate the origin of 
rotons in a Bose system. 



In physics the most fascinating phenomena are perhaps 
so-called macroscopic quantum phenomena, such as su- 
perfluidity, partially due to that they are anti-intuitive 
at a visible scale. It might look unHkely to get a clear 
theoretical picture of these phenomena since that quan- 
tum many-body problems are complex problems. How- 
ever, some direct analysis still can be performed and some 
understanding can be obtained at the very fundamental 
level, for example, the construction of Laughlin's wave- 
function for fractional quantum Hall effect. In this paper, 
we attempt to argue that supersolidity [l| is a natural 
consequence of Bose exchange symmetry through a di- 
rect analysis of the many-body wave functions of a Bose 
system, our analysis also reveals the origin of rotons in 
a Bose system, which is another fascinating issue in con- 
densed matter physics. 

The knowledge of superfluidity has accumulated 
through the efforts of several generations of physicists 
Landau first related superfiuidity with the prop- 
erties of the spectrum of a system, specifically, the quasi- 
particle spectrum. Bloch's Q and Leggett's works 
clearly point out that a more transparent way to un- 
derstand this phenomenon is based on the properties of 
many-body dispersion spectrum, i.e., the lowest eigen en- 
ergies of the system at given momenta. When there exist 
local minima in the dispersion spectrum, the energy bar- 
riers which separate the minima will prevent the decay of 
current carried by the states corresponding the minima 
(see Fig. [2] for example), thus superfiuidity occurs . 

We have argued in [a| that the existence of the local 
minima in the dispersion spectrum of a Bose system could 
be understood in terms of Bose exchange symmetry. Let 
us consider a 2D periodic Bose system with N particle, 
with a finite-range repulsive interaction and with a square 
geometry (see Fig.[T]), The Hamiltonian takes a form 
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Where M is the mass of a particle, g is the interaction 
strength, / describes the form of a repulsive interaction, 
2ttR is the linear size of the system. 9i = {xi,yi)/R 
where {xi,yi) are the coordinates of the z-th particle. We 
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Figure 1: The geometry of a 2D periodic system. The edge 
AB is identified with DC and BC identified with AD. The 
linear size of the system is 2-nR. 
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Figure 2: The spectra of yrast states with Py = 0, relative to 
the ground state energy, for three systems with N = 9, and 
having various interaction strengths g. the interaction takes 



a form of /(6») = An^ 



/(47r^^) with C = 27r/15. 
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will adopt the unit system generated by h, R and M, 
except that the energy is in the unit of tf/2MR'^ where 
a factor of a half is involved for convenience. 

The many-body spectrum can be labeled by the mo- 
mentum of the system P = {Px,Py), since Px,Py are 
good quantum numbers due to the translational symme- 
try. We consider part of the spectrum, E = E{Px,Py), 
where E is the minimum energy at the given {Px, Py). We 
will call the many-body states corresponding this disper- 
sion spectrum the yrast states, following convention in 
nuclear physics. 

Due to Galileo invariance [1, B] , we can derive the full 
yrast spectrum if part of this spectrum in the regime 
< Px < N, < Py < N is known. Particularly, we 
can derive all local minima of the spectrum and realize 
that there is an upper momentum limit beyond which no 
supercurrents exist any more, i.e., there exists a critical 
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velocity |J, [6| . We will only focus on the structure of the 
yrast spectrum in the momentum regime < < iV, 
< Py < N in the rest of this paper. 

The perturbation analysis of the yrast states 
and dispersion spectrum of the system at small 
g, is straightforward and presented in For 
example, the yrast states for < Px < N, 
Py = are be approximated well by Fock states 
1(0,0)^(1,0)0), 1(0,0)^-1(1,0)1), 1(0,0)^-2(1,0)2), 

1(0,0)0(1,0)^), where {m,n), a pair of inte- 
gers, denotes the single particle orbit ^m,n{x,y) = 
l/(27r)e™^/'"+*"2'/-« and |(0, 0)^i (1, 0)^^) denotes a Fock 
state with ji particles occupying orbits (0, 0) and j2 
particles in orbits (1,0). These Fock states will have 
the same (direct) Fock energy while the exchange 
energy, as a function of Pa;, has a paraboHc behavior, 
Eex oc g{N — Px)Px- This behavior of the exchange 
energy, without any classical analog, is responsible for 
the possible existence of local minimum of the dispersion 
spectrum at {Px,Py) = (-/V, 0). In we also present 
the second microsmany-body argument for the exchange 
origin of superfluidity. By comparing the dispersion 
spectrum of a spinless Bose system and that of a 
corresponding spinor Bose system, one can naturally 
imagine that the possible energy barriers separating two 
neighboring minima of dispersion spectrum (of a spinless 
system) will become more apparent with increasing g, 
i.e., strong interaction enhances superfluidity. 

In this paper, we will present another many-body ar- 
gument for the quantum origin of superfluidity, which il- 
lustrate in a direct way the role of exchange symmetry in 
determining the dispersion spectrum of the systems. This 
argument also indicate a natural compatibility of super- 
fluidity and crystaUine ordering, on the contrast to previ- 
ous common beliefs that crystalline ordering will destroy 
superfluidity. Supersolidity is a natural consequence of 
this compatibility. 

We consider the system specifled above with = 9, 
the generalization for N = m? , where m is an integer 
is straightforward. The dispersion spectra shows that 
with large g there are also supercurrents at {Px,Py) = 
(3,0), (6,0) besides the one at {Px,Py) = (9,0) (see 
Fig. [2]). The supercurrent states at P = (3,0) and (6,0) 
are rather surprising, in which it is impossible that all 
nine particles occupy the same orbits. This 'new' super- 
currents can be explained. 

A many-body state of the system can be written in the 
spatial presentation 



E 

xi,X2,---,xa 



tp{xi,X2, ...,X9)\Xi,X2, 



,X9) 



(2) 



A1 


A2 


A3 


X 


X 


X 


A4 


A5 


A6 


X 


X 


X 


A7 


A8 


A9 


X 


X 


X 



A1 


A2 


A3 


X 


X 


X 


A4 


A5 


A6 


X 


X 


X 


A7 


A8 


A9 


X 


X 


X 



A1 


A2 


A3 


X 


X 


X 


A4 


A5 


A6 


X 


X 


X 


A7 


AS 


A9 


X 


X 


X 



Figure 3: Some examples of configurations with 27r/3 — x — 
symmetry and with rather small interaction energies. Crosses 
mark the positions denoted by AiS. 



and the interaction energy of the state has the form of, 

N 



one might refer to each state \xi,X2, ■■■,xq) as a config- 
uration with the first particle at position Xi, the second 
at X2, the ninth at xg. 



eL = mY.f{x,-x,m (3) 

N 

= E H{xl,...,Xg)\^g{Y,f{x^-' Xj)) (4) 

Xi,...,Xg i<j 

which can be interpreted as average of the interaction en- 
ergy of each configuration, with weight being the square 
module of the wavefunction. The configurations with 
small interaction energy shall be largely explored by the 
yrast states at large g, in which the interaction energy 
dominates the kinetic energy. 

We shall consider a group of configurations with a cer- 
tain discrete translational symmetry and its relationship 
with yrast states. The symmetry is following. If one 
configuration is translated by a mount of 27r/3P along 
X axis, the first particles will be in the previous posi- 
tion of the second particle {xi X2), the second par- 
ticles will reach the previous position of the third par- 
ticle {x2 X3), and X3 ^ Xi, Xi —>■ x^ ^ xe ^ x^, 
X7 ^ xs ^ Xf) ^ xi. We label the positions of 
nine particles of this configuration hy Ai, A2, Ag, i.e., 
xi = Ai,X2 — A2,...,xg = Ag. Thus after the transla- 
tion, Xi = A2,X2 = A3, X3 = Ai, ...,Xg = At. 

We refer to this symmetry as 27r/3 ~ x — symmetry. 
Obviously, most configurations with 2-^-/3— x— symmetry 
has a relatively small interaction energy since the parti- 
cles are well separated along x axis, (see Fig.[3]for exam- 
ples). 

We shall argue that the configurations with 27r/3 — a: — 
symmetry can only be accommodated (with a non-zero 
weight) by yrast states at Px — 0, 3, 6, 9, ... . 

First, with this translation, the wavefunction associ- 
ated with the configurations changes to be, 

«'(X1 = A2,...,Xg = A7) = i;{Ai+b,...,Ag + b) 

= e*''^/V(^i,A2,...,^) (5) 

Where b = (27r/3P, 0), P — {Px,Py) is the momentum 
operator. For eigenstates of the P, bP/h ~ Pr/S * 27r 
with Px in the unit we specified above. 

Due to Bose exchange symmetry, the wavefunction is 
invariant under any permutation of the coordinates of 
particles, thus 
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V'(A2,A3,...,^7) =^(^l,^2,-,^). (6) 

Combining Eqs.El E one is led to 





l)ip{Ai,A2,A3;A4„A5,A6 



A7,As,Ag) = 0. 

(7) 

Thus ip{Ai, A2, A3, Ai, A5, Ae, At, As, Aq) ^OatP^^ 
0,3,6,... . The configurations with 27r/3 — a; — symmeirj/ 
can't be effectively involved in those yrast states with 
P. ^0,3,6,... . 

One can also note that if one configuration is not in- 
volved in a yrast state jV"), its neighboring configura- 
tions can't be effectively involved in \ip). For a con- 
figuration with first particle at position Bi, the second 
particle in position B2, ■■■ , in its neighboring config- 
urations, the first particle is at position Bi + Si where 
Si is small position shift, and the second particle is at 
-82 +(^2 with (52 being small, ... . li 'ip{Bi, B2, Bg) — 0, 
then lipiBi + 5i, B2 + S2, Bg + Sg)\ can't be large, oth- 
erwise the value of the wavefunction changes dramat- 
ically from {xi,X2,---) = i?2, •••); to {xi,X2,---) = 
{Bi+Si,B2 + S2, ...), which means that the wavefunction 
has a large gradient and thus a large kinetic energy (den- 
sity) near this neighboring configuration, but large ki- 
netic energy (density) shall be avoided in the yrast states 
which pursues the minimum energy at given momenta. 

Therefore, the yrast states at Px — 0, 3, 6, 9, ... can use 
the configurations with 2^/3 — x — symmetry and their 
neighboring configurations to lower their interaction en- 
ergy while the yrast states at other Px can not. Similarly, 
the yrast states at Py = 0, 3, 6, 9, ... can effectively access 
configurations with 2Tr/3—y— symmetry and their neigh- 
boring configurations to lower the interaction energy 
while the states at other Py can not. At sufficient large g, 
the yrast states at {Px.Py) = (3, 0), (6, 0), (9, 0), (3, 3), 
will become supercurrent states, with energies lower than 
the energies of their neighboring yrast states. 

The configurations with 27r/3 — x,y — symmetry are 
also the fundamental components to build an inner 
crystal-like structure in the system, one thus draw a con- 
clusion that only those yrast states with particular mo- 
mentum values can support such inner crystal-like struc- 
tures, they are supercurrent states. The inner crystal-like 
structure of the ground state at large g can be seen clearly 
in the plot of pair correlation function (see Fig.|4l), which 
is defined as 



where 6a is the referees point. 

The above discussions can be generalized to systems 
with large N and to 3D systems straightforwardly. In 
a 2D system with square periodic boundary geometry, 
with A'' = Q, and with large g, the yrast spectrum 



Figure 4: Pair correlation profile of the ground state of a 
system with 5=1 (left) and of a system with g = 12 (right). 
The reference point 0a is located at the center. N = 9. At 
rather small g, g = 1, the system has short-range correlations 
but no long-range correlations and is in a gas- like or liquid-like 
phase (left). At large g, g ~ 12, the system has long-range 
correlations and is in a solid- like phase (right). 
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Figure 5: The yrast spectrum of a 3D periodic system, rel- 
ative to the ground state energy. Py — Pz — Q. N — 8. the 

interaction is V(0) = (27r)^ffe~*'''''«V(V4^0^ with g = 6 , 
^ = 37r/20, where = {9x, dy,dz) is the difference between the 
generalized angular coordinates of two interacting particles. 



will has local minima at {Px,Py) — (m, 0), (2to, 0),..., 
(m,m), {2m,m),...,{m,m), due to the compatibility 
of the yrast states at these momenta with the particle 
configurations with 27r/m — x,y — symmetry. In a 
3-D system with cubic periodic boundary geometry 
and with A^ = m^ , the yrast states at {Px, Py, Pz) = 
(m, 0, 0), {2m, 0, 0),..., (m, to, 0), (2m, m, 0), ...,{m, to, to), 
{m^,m^,m^) can become supercurrent states 
for their accommodations of the configurations with 
21: jm — x,y,z — symmetry. Numerical results of a 3D 
periodic cubic system with A^ = 8 is plotted in Fig. [5] . 

The above argument for exchange origin of superfiuid- 
ity can also be applied to a system with an interparticle 
interaction which is repulsive at short range and attrac- 
tive at long range, for example, the Helium system. In 
these systems, the configurations with that type of dis- 
crete translation symmetries are still energetically favor- 
able and the states compatible with these configurations 
can become supercurrent states. One might naturally 
speculate that supersolidity can be observed in neon sys- 
tems, much like the case of helium systems Q. 

It is interesting to note that the smallest momentum of 
a supercurrent in a Bose system with large repulsive in- 
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teraction is (5p = N^l'^ (d is the dimensionality) and that 
it corresponds to the inverse of the average particle dis- 
tance a. In liquid ^He, the momentum corresponding to 
rotons in quasiparticle spectrum is also around hja (h is 
Planck's constant), which naturally suggests that rotons 
are due to the local minimum of the yrast spectrum. The 
quasiparticle dispersion curve lies above the many-body 
dispersion curve, as illustrated schematically in Fig. [6l 
Corresponding to the local minima of the yrast spectrum 
at momentum hja, there are rotons in the quasiparticle 
spectrum. 



h/a 



Figure 6: Schematic plots of the quasiparticle dispersion (solid 
line) and yrast spectrum (dashed line) of a repulsive Bose sys- 
tem. Quasiparticles generally do not exist beyond a thresh- 
old where certain decay channels become allowed h is 
Planck's constant. 



In conclusion, we present a microscopic many-body ar- 
gument for supersolidity. We also predict possible super- 
solidity in ^"Ne and in ^^Ne. 



[1] E. Kim and M. H. W. Chan, Nature 427, 225 (2004); E. 
Kim and M. H. W. Chan, Science 305, 1941 (2004);E. 
Kim and M. H. W. Chan, Phys. Rev. Lett. 97, 115302 
(2006); E. Kim and M. H. W. Chan, Low. Temp. Phys. 
138 (3-4), 859 (2005); A. C. Clark, J. T. West and M. 
H. W. Chan, Phys. Rev. Lett. 99, 135302 (2007); X. Lin, 
A. Clark and M. H. W. Chan. Nature 449, 1025 (2007). 
For previous theoretical predictions, see M. Wolfke, Ann. 
Acad. Sci. Techn, Varsovie 6, 14 (1939); W. H. Keeson, 
Helium, p. 374, Elsevier, New York, 1942; A. F. Andreev 
and I. M. Lifshitz, Sov. Phys. JETP 29, 1107 (1969); G. 
V. Chester, Phys. Rev. A 2, 256 (1970); A. J. Leggett, 
Phys. Rev. Lett. 25, 1543 (1970). 

[2] see D. R. Tilley and J. Tilley, superfluidity and supercon- 
ductivity, and referees therein, 3rd ed., Adam Hilger Ltd, 
Bristol, 1990. 

[3] L. D. Landau, J. Phys. USSR 5, 71 (1941). 

[4] F. Bloch, Phys. Rev. A 7, 2187 (1973). 

[5] A. J. Leggett, Rev. Mod. Phys. 73, 307 (2001). 

[6] Y. Yu, cond-mat/0609711 

[7] One notes that the system can form an energetically fa- 
vorable crystal only at A'' = m^. If the system has differ- 
ent periodic boundary geometry and/or different sizes, it 
will form a crystal at different N values. Since a crys- 
tal structure is expected at strong interaction regime, we 
thus consider only the systems with these specific A'' val- 
ues. 

[8] Here one ignores the probability for a supercurrent state, 
a yrast state corresponding to a local minimum, to decay 
into another supercurrent state (with a less energy). In 



such a decay process, the system emits a quantum with 
a certain energy and with a large momentum which is 
at least of the order of h/a where a is the average inter- 
particle distance. Depending on its own low energy spec- 
trum, the environment may not be capable of absorbing 
this quantum and thus this decay process is unlikely to 
occur. 

[9] Supersolidity was searched in hydrogen systems but was 
not observed, see A. C. Clark, X. Lin, and M. H. W. 
Chan, Phys. Rev. Lett. 97, 245301 (2007). A solid hydro- 
gen, however, might not be considered as an exact spin- 
less Bose system due to the conversion between ortho-Il2 
and para-H2 in it. 

[10] In a 2D or 3-D periodic Bose system, both the yrast spec- 
trum and the quasiparticle dispersion are anisotropic. If 
one ignores the anisotropy of the quasiparticle dispersion 
and considers the average of the lowest quasiparticle en- 
ergies at P = |P|, then this average is expected to have 
a minimum around the momentum h/a. 

[11] L. P. Pitaevskii, Sov. Phys. JETP 9, 830 (1959); L. 
D. Landau and E. M. Lifshitz Courses of Theoretical 
Physics, Vol. 9, Statistical Physics Part 2 (by E. M. 
Lifshitz and L. P. Pitaevskii) 125-139 (Pergamon, Ox- 
ford, 1981); A. A. Abrikosov, L. P. Gorkov and I. E. 
Dzyaloshinski, Methods of quantum Field theory in sta- 
tistical Physics 237-247 (Dover Publications, New York, 
1975); M. B. Stone, I. A. Zaliznyak, T. Hong, C. L. Bro- 
holm and D. H. Reich, Nature 440, 187 (2006). 



